Visualizing the Environmental Impact of Beef Consumption using Plotly and Shiny

##constants
water.rate <- 6.36 #(million gallons of water per 1 animal) 
co2.emision <- 102.96 #(thousand lbs of co2 emissions per 1 animal) 
land.use <- 77 #(acres of land per 1 animal)

##load in the dplyr library
library(dplyr)

##entry weight for pound of meat
weight = 100

## calculations for enviornmental impact  
carcass_calculator_data_100_pounds <- carcass_calculator_data %>%
  mutate(number_cows = ceiling(weight/total_weight)) %>% 
  select(-total_weight)

##calculate the enviornmental impact 
carcass_calculator_data_100_pounds <- carcass_calculator_data_100_pounds %>% 
  mutate(water_rate = water.rate * number_cows) %>% 
  mutate(co2_emission = co2.emision * number_cows) %>% 
  mutate(land_use = land.use * number_cows) 

## load in the required the packages 
library(ggplot2)
library(wesanderson)

##create a color palette
values=wes_palette(n=4, name="Darjeeling1")

##reorder the data in order of number of cows 
plot.data <- carcass_calculator_data_100_pounds %>% 
  mutate(cut = reorder(cut, number_cows))

##plot th data in a barplot 
plot.out <- ggplot(plot.data, 
                   aes(x = cut, 
                       y = number_cows)) + 
  geom_bar(stat="identity", fill = values[4]) + 
  coord_flip() + 
  xlab('') + 
  ylab('Animals') + 
  theme(legend.position="none",
        axis.title=element_text(size=12),
        plot.title=element_text(size = 16)) + 
  ggtitle(paste0('Animals for ', weight, ' Pounds of Meat by Cut'))

## make a vector of a subset of the cuts  
cut.list <- c('Tongue',
              'Ground Beef',
              'Tenderloin',
              'Bone-in Ribeye',
              'Skirt')

## plot the data for the subset of the cuts
plot.out %+% filter(plot.data, cut %in% cut.list) + 
  theme(axis.text=element_text(size=14))

library(plotly)
ggplotly(plot.out, tooltip = c("x", "y"))

x <- list(
  title = "Distance (in miles)"
)

y <- list(
  title = "Elevation (in feet)"
)

plot_ly(data = devils.path,
        x = ~Distance..miles.,
        y = ~Elevation..feet.,
        type = 'scatter', 
        mode = 'lines',
        color = I('purple'),
        hoverinfo = 'text',
        text = paste(
          round(devils.path$Distance..miles., 1), 
          "miles, ",                    
          round(devils.path$Elevation..feet., 0),
          "feet elevation"),
        fill = 'tozeroy')  %>%
  layout(xaxis = x, 
         yaxis = y,
         title = "The Devil's Path",
         shapes=list(type='line', 
                     y0= 3500, 
                     y1= 3500, 
                     x0=min(devils.path$Distance..miles.),
                     x1=max(devils.path$Distance..miles.),
                     line=list(dash='dot', width=1)))